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Abstract: 

We have studied the effect of one-loop logarithmic electroweak radiative corrections 
on WZ and W7 production processes at the LHC. We present analytical results for the 
leading-logarithmic electroweak corrections to the corresponding partonic processes du — 
WZ,W7. Using the leading-pole approximation we implement these corrections into 
Monte Carlo programs for pp — ^ luil'l', Ivi'y. We find that electroweak corrections lower the 
predictions by 5-20% in the physically interesting region of large transverse momentum 
and small rapidity separation of the gauge bosons. 
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1 Introduction 



Vector-boson pair production provides us with an important testing ground for the non- 
abehan structure of the Standard Model (SM). While gauge-boson properties, such as 
masses and couplings to fermions, have already been measured with great accuracy at 
LEP and Tevatron, vector-boson self-interactions have not been tested with comparable 
precision. New physics occurring at energy scales much larger than those probed directly 
at forthcoming experiments could modify the structure of these interactions. These mod- 
ifications are parametrized in terms of anomalous couplings in the Yang-Mills vertices. 

In the last few years, the contribution of trilinear gauge-boson couplings was directly 
measured via vector-boson pair production at LEP2 and Tevatron. While, in particular, 
LEP2 has been able to produce W"'"W~ pairs with high statistics, nevertheless all these 
events were generated at rather modest centre-of-mass (CM) energies (-Ecm ~ 210 GeV). 
The effect of anomalous couplings is, on the other hand, expected to be strongly enhanced 
by increasing the invariant mass of the gauge-boson pair Myv' iY^y = W, Z,7), since 
these couplings in general spoil the unitarity cancellations for longitudinal gauge bosons. 
Hence, at future colliders it will be useful to analyse the di-boson production at the highest 
possible CM energies. 

Moreover, vector-boson pairs constitute a background to other kinds of new-physics 
searches. One of the gold-plated signals for supersymmetry at hadron colliders is chargino- 
neutralino pair production, which would give rise to final states with three charged leptons 
and missing transverse momentum [[![]; the primary background to this signature is given 
by WZ or W7* production. Leptonic final states, coming from WZ or W7*, could also 
fake WZ vector-boson scattering signals, or W^W"^ and W^W^ scattering signals if one 
of the charged leptons is lost in the beam pipe, which are again expected to be enhanced 
at high CM energies [ ^ . 

In the next years, hadron colliders will be the main source of vector-boson pairs with 
large invariant masses Myy'- Tevatron Run II will collect from tens to hundreds of 
events, depending on the particular process. The large Hadron Collider (LHC) will further 
increase the event number by roughly two orders of magnitude [ ^ . Owing to the expected 
increase in statistics, theoretical predictions must reach high accuracy to allow for a decent 
analysis of the data. 

Hadronic di-boson production has received a lot of attention (for a review on the sub- 
ject see Ref. [ ^). Originally computed by treating W and Z bosons as stable particles, 
tree-level cross sections for W'^'W", W^Z, ZZ, W^7, and Z7 production and decay have 
been updated by evaluating, in narrow-width approximation but retaining spin informa- 
tion via decay-angle correlations, the doubly-resonant contribution to the four-fermion 
final states (e.g. qq' W^Z 4/) [ Q and the resonant contribution to two-fermion 
plus photon final states (e.g. qq' W^7 — 2/ + 7) [|^. As a further step, in the last few 
years Monte Carlo programs [ Q have included the full qq' —>■ 4/ amplitude, by taking into 
account finite-width effects and the irreducible background owing to non-doubly-resonant 
diagrams. 

The 0{as) QCD corrections to gauge-boson pair-production and decay have been ex- 
tensively analysed by many authors. Gauge-boson pair-production cross sections have 
been calculated at next-to-leading order (NLO) accuracy retaining the full spin corre- 
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lations of the leptonic decay products. Several NLO Monte Carlo programs have been 
implemented and cross checked so that complete 0{as) corrections are now available [ 
H, 1^, p. QCD corrections turn out to be quite significant at LHC energies. They can 
increase the lowest-order cross section by a factor two if no cuts are applied and by one 
order of magnitude for large transverse momentum or large invariant mass of the vector 
bosons [ 1^, §. By including a jet veto, their effects can be drastically reduced to the order 
of tens of per cent [ ^, but in any case they have to be considered to get realistic and 
reliable estimates of total cross sections and distributions. 

In view of the envisaged precision of a few per cent at the LHC, also a discussion 
of electroweak corrections is in order. For single W- and Z-boson production, 0{a) 
corrections have been computed taking into account the full QED and weak contributions 



[ |T0[- For gauge-boson pair production at hadron colliders, the electroweak corrections 
have been taken into account only via an effective mixing angle. 

As well known, the impact of 0{a) electroweak contributions grows with increasing 
energy. Analyses of the high-energy behaviour of electroweak corrections in general and 
for specific e+e~ and 77 processes have already been performed revealing effects which 



should be clearly visible at future linear colliders (see for instance Refs. [ ^ Q)- At high 



energies the electroweak corrections are dominated by double and single logarithms of the 
ratio of the energy to the electroweak scale. In Refs. [ Q it has been shown that 



the leading-logarithmic one-loop corrections to arbitrary electroweak processes factorize 
into the tree-level amplitudes times universal correction factors. These results represent 
a process-independent recipe for the calculation of leading logarithmic corrections. 

Using the method of Refs. [13, n|, we investigate in this paper the effect of leading- 



logarithmic electroweak corrections to the hadronic production of W^Z and W^7 pairs 
in the large-invariant-mass region of the hard process at the LHC. Since the aim of this 
paper is to describe the structure of the 0{a) electroweak corrections and to give an 
estimate of their size, we have not included QCD corrections. Also, QED corrections 
are not fully considered as they strictly depend on the experimental setup. We omit all 
infrared- singular terms originating from the massless photon, as explained in App. |A.l 



and focus on the contributions of the leading electroweak logarithms originating from 
above the electroweak scale. 

The simplest experimental analyses of gauge-boson pair production will rely on purely 
leptonic final states. Semi-leptonic channels, where one of the vector bosons decays 
hadronically, have been analysed at the Tevatron [ ^ showing that these events suf- 
fer from the background due to the production of one vector boson plus jets via gluon 
exchange. For this reason, we choose to analyse only di-boson production where both 
gauge bosons decay leptonically into e or fi. 

The paper is organized as follows: in Sect. § we give some details on the general setup 
of our calculation. In Sect. |^ and App. ^the logarithmic electroweak one-loop corrections 
are examined and presented in analytical form. Section ^ contains a numerical discussion 
for WZ production and decay, while Sect, ^covers W7 production and decay. Our findings 
are summarized in Sect. ^. 
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2 Processes and their computation 



We consider in detail two classes of processes, 

(i) pp — > luil'l', with /, /' = e, /i, and 

(ii) pp /z/;7, with I = e, fi. 

The first class is characterized by three isolated charged leptons plus missing energy in 
the final state. In our notation, lui indicates both and This kind of processes 
includes WZ production as intermediate state. The second class is instead related to W7 
production. Both classes of processes are described by the formula 

da'^^'^^(Pi,P2,P/) = E/ da;idx2 (xi, Q2)/,,/,,(x2, Q') da*^(xiPi, xsPs,^/), (2.1) 

where pf summarizes the final-state momenta, fi^hi and fj^h2 are the distribution functions 
of the partons i and j in the incoming hadrons hi and /i2 with momenta Pi and P2, 
respectively, Q is the factorization scale, and a*-' represent the cross sections for the 
partonic processes. Since the two incoming hadrons are protons and we sum over final 
states with opposite charges, we find 

da"^"^(Pi,P2,p/) = Jdxidx2 E E [/D,p(a;i,Q')/u,p(x2,Q')da°U(a;iPi,x2P2,P/) 

U =u,c D=d,s 

+ /u,p(a;i,<5^)/D,p(a;2,Q^) da^°(xiPi,X2P2,P/) 
+ /D,p(a;i, Q')/u,p(a;2, Q') da°U(a;,Pi, x2P2,Pf) 
+ fvA^u Q')hA^2, Q') daUO(a;iPi, X2P2,Pf)] (2.2) 

in leading order of QCD. 

The tree-level amplitudes for the partonic processes have been generated by means of 



PHACT [ |T^, a set of routines based on the helicity-amplitude formalism of Ref. [|T7 . 
For the numerical results presented here, we have used the fixed-width scheme with 
Tz = 2.512 GeV and Tw = 2.105 GeV, and the input masses Mz = 91.187 GeV and 
M-w = 80.45 GeV. The weak mixing angle is fixed by = 1 — M-^/M|. Moreover, 
we adopted the so called G^-scheme, which effectively includes higher-order contributions 
associated with the running of the electromagnetic coupling and the leading universal 
two-loop mt-dependent corrections. This corresponds to parametrize the lowest-order 
matrix element in terms of the effective coupling ac^ = v^G^jM-^s^/tt. However, we 
use a(0) = 1/137.036 for the coupling of the real photon in pp li^ij, i-e. we multi- 
ply the corresponding cross sections in the G^-scheme by a{0)/ac^- Additional input 
parameters are the quark-mixing matrix elements whose values have been taken to be 
\Vud\ = |Ks| = 0.975, |Vus| = iKdl = 0.222, and zero for all other relevant matrix ele- 
ments. 

As to parton distributions, we have used CTEQ(5M1) [ [18| at the factorization scales 
Q2 = i (m^ + Mi + P^ilui) + P^l'l')) (2.3) 
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and 

Q' = l{M^ + P^{lui) + PU^)) (2.4) 

for WZ and W7 production processes, respectively, where Pt denotes the transverse 
momentum. This scale choice appears to be appropriate for the calculation of differential 
cross sections, in particular for vector-boson transverse-momentum distributions [ |§, ^ • 

We have, moreover, implemented a general set of cuts, proper for LHC analyses, 
defined as follows: 

• lepton transverse momentum Pt(0 > 20 GeV, 

• missing transverse momentum Pr^^^^ > 20 (50) GeV for WZ (W7), 

• lepton pseudo-rapidity \rii\ < 3 (2.5) for WZ (W7), where rji = — log (tan^^.j/2), 9i is 
the polar angle of particle i with respect to the beam, and i = 

• rapidity-azimuthal-angle separation ARi^ = ^ {rji — rj^y + {(pi — (p^y > 0.7 be- 
tween charged lepton and photon for W7. 

For the different processes considered, we have also used further cuts which are described 
in due time. In the following sections, we present results for the LHC at CM energy 
y/s = 14TeV and an integrated luminosity L = 100 fb~^. 



3 Electroweak 0(cx) corrections 

We are interested in the electroweak 0{a) corrections to the processes pp Ivil'V and 
pp lui'j in the region of phase space where these are dominated by the gauge-boson 
pair-production subprocesses pp — > WZ and pp W7, respectively. In this region, the 
dominant contributions are those that are enhanced by the resonant propagators of the W 
boson and in the first process also of the Z boson. These can be most effectively calculated 
in the so-called leading-pole approximation (LPA), which is a double-pole approximation 
(DPA) for pp WZ — s> luil'l' and a single-pole approximation (SPA) for pp — > W7 — > 
li'i'j. The LPA has been successfully applied for the calculation of electroweak corrections 
to W-pair production [ |TP|, ^, gg. 

At tree level, the DPA for the partonic process qq' — > WZ — luil'l' reads 



M 



Born,DPA p^, - + iMwFw p| - M| + iMzFz 

X E Mli-'^^'-MZ^r^-Ml^-'', (3.1) 

A,A' 

where the gauge- dependent doubly-resonant contribution is replaced by the well-defined 

gauge-independent residue, and A^Bom^^^^'i -^BotiT''^'' ^^nd A^Born*' ' denote the on-shell 
Born matrix elements for the boson production and decay processes. The sum runs over 
the physical heli cities A, A' = 0, ±1 of the on-shell projected W and Z bosons (see App. A 
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of Ref. [ 2^ for details), while the momenta of the virtual W and Z bosons are denoted 
by pw and pz, respectively. For the process qq' W7 — * /z//7, the SPA is given by 

In LPA, the 0{a) electroweak corrections to boson-production processes can be di- 
vided into two classes, factorizable and non-factorizable corrections. The non-factorizable 
corrections, i.e. those contributions that cannot be associated with either boson pro- 
duction or decay, have been evaluated for boson-pair production in e+e" annihilation 
in Refs. [ |2^, There, these corrections turned out to be small. We assume that this 
holds as well for the similar processes considered here and do not consider non-factorizable 
corrections any further. 

Moreover, we do not include real corrections and restrict our discussion to the infrared- 
finite part of the virtual factorizable corrections, as defined in App. |A.1| . These contri- 
butions can be expressed in terms of the corrections to the on-shell boson production 
and decay subprocesses. The matrix element for the virtual corrections to the process 
qq' —>■ WZ —y Ivil'l' can be written as 



^-^virt.DPA 



pIj-M^ + iMwFw vl - Ml + iMzFz 

2^1 virt -'^'Born •'^'Born 

A,A' 

+ All-:^^^^'A<^o^-'^'5A<.^r''''}, (3.3) 

where (5A^^^j,t~^^^^^', 5A^^X^^^\ and dM^i-^.^^ ' denote the virtual corrections to the on- 
shell matrix elements for the boson production and decay processes. A similar expression 
holds for the process qq' — >• W7 Ih'i'j. 

We focus in particular on the corrections involving single and double enhanced 
electroweak logarithms at high energies, i.e. on 0{a) contributions proportional to 
Q;log^(s/M-^) or alog(s/M-^), where \/S is the CM energy of the partonic subprocess. 
The logarithmic approximation yields the dominant corrections at CM energies large com- 
pared to the gauge-boson masses, s S> M^. In the high-energy limit, however, there might 
be also enhanced non- logarithmic contributions that are a priori relevant. In general they 
contain constant terms proportional to M^/M^ and m^/M^. While for transverse gauge 
bosons, these Higgs- and top-mass-dependent corrections are entirely due to renormaliza- 
tion effects and can be effectively accounted for by using the G^-scheme, for longitudinal 
gauge bosons additional contributions of this kind exist. These non-logarithmic 0{a) con- 
tributions are process-dependent. For e"''e~ W^W~, where complete 0{a) corrections 
and their high-energy limit are available [ |ll|, the above terms turn out to be of order of 



a few per cent. We can then qualitatively assume that this holds as well for similar pro- 
cesses like hadronic di-boson production, even if only an exact computation could really 
furnish a precise statement on this point. Neglecting non- logarithmic terms can therefore 
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be considered a reasonable approximation at the LHC, where the experimental accuracy 
in the high-energy regime is at the few-per-cent level. 

Since the decay processes involve no large-energy variable, the corresponding virtual 
corrections vanish in the logarithmic approximation. As a consequence, we do not consider 
in the following the last two contributions on the right-hand side of (|3.3| ). Moreover, for 
the boson production processes qq' —>■ WZ, W7 we take into account only the correction 
to the dominating channels involving two transverse (TT) or two longitudinal (LL) gauge 
bosons. The contributions of the mixed (LT, TL) channels are suppressed relatively 
to the others by factors of Mw/\/^ in the high-energy limit (see Fig. and thus the 
corresponding corrections are unimportant. 

The logarithmic virtual electroweak corrections to the dominating channels are cal- 
culated using the general method given in Refs. [O, KM. The corresponding analytical 



expressions for the processes du — >■ W^A^a, -A = Z,7, are given in App. ^ Those for 
du —>■ W^Nx are derived via CP symmetry [see ( [A. 21 )]. Our predictions are obtained by 
considering the matrix element squared 



\M\''= |A^Bornl' + 2Re 



(3.4) 



where A^Bom is the exact Born amphtude, while the 0{a) contribution is computed in 
LPA based on ( p.l|) - (|37^ ), and the formulas given in App. ^ In the high-energy limit, a 
reasonable approach is to neglect fermion and boson masses, as compared with V^, wher- 
ever possible. The expressions given in App. ^ are based on this approximation. However, 
we take into account the exact kinematics by evaluating the complete four-fermion or two- 
fermion-plus-photon phase space and use the exact values of the kinematical invariants 
in all formulas. Moreover, we do not use the high-energy approximations ( A. 7 ) and (|A.8|) 



in the correction factors but we implement the 0{a) contributions according to the full 
expressions given in (|Al2|) , ( [AlsD , ([AlsD , and ( |^:21| )- (|A:23|) with the exact (SU(2)- 



transformed) Born matrix elements. Owing to our choice of the input-parameter scheme, 
the terms proportional to Aa(M^) in ( |A.22| ) and ( |A.23| ) are omitted, since these are 
already taken into account by using ac^ instead of a{0) as input. 

In the universal logarithmic corrections given in App. ^ the pure angular-dependent 
logarithms, such as a log^(|f|/s) and alog(|f|/s) with f equal to the Mandelstam variables 
t and u of the partonic subprocess, are not included. The validity of this approximation 
relies therefore on the assumption that all the variables s, \t\, and l^l are large compared 
with and approximately of the same size, 

s ~ |t| ~ |m| > Mw. (3.5) 

This implies that the produced gauge bosons have to be emitted at sufficiently large an- 
gles with respect to the beam. Hence, the validity range of the high-energy logarithmic 
approximation for the radiative corrections corresponds to the central region of the boson 
scattering angle in the di-boson rest-frame. For s-channel processes, integrating over the 
full angular domain does not affect the reliability of the result at logarithmic level, since 
the neglected pure angular- dependent logarithms would give rise only to subleading con- 
stant terms, if included. For t-channel dominated scatterings like WZ or W7 production, 
the situation is instead more delicate. The t-channel pole in the Born matrix element 
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gives rise to additional enhanced logarithms when integrated over the full kinematical 
range. Since these terms are not included in our 0{a) analysis, we have to take care 
that we do not get sizeable contributions from small scattering angles with respect to the 
beam. On the other hand, our formulas do not fake spurious contributions as long as 
s, |u| ^ M^, since the large logarithms become small for s, |t|, \u\ ~ M^. 

4 W^Z production 

In this section, we present some cross sections and distributions for the leptonic processes 
pp —>■ Ivil'V with = e, /i. These final states allow to analyse WZ production and thus 
in particular to test trilinear gauge-boson couplings. Systematic studies of the effect of 
anomalous couplings on the hadronic production of gauge-boson pairs have shown that 
deviations from the SM cross sections should be particularly enhanced when gauge bosons 
are produced at high CM energies and at large scattering angles in the di-boson rest-frame. 
The same kinematical region is also proper to search for scatterings of strongly interacting 
vector bosons. 

It is therefore particularly interesting to study the electroweak corrections in these 
kinematical configurations, where their effect is also expected to be more sizeable. As 
an illustration of the behaviour and the size of the 0{a) contributions, we have chosen 
to analyse the distribution of the reconstructed Z-boson transverse momentum Pt{1'1'). 
The Pt variable is commonly used at hadron colliders because large Pt requires high 
CM energies and large angles. We study also pure angular observables of interest in the 
high-energy regime of the hard scattering. 

4.1 Born level 

We start recalling basic properties of the Born amplitude, which are useful later in dis- 
cussing radiative-correction effects. In Fig. |l], just for explicative purposes, we have plot- 
ted the on-shell Born cross section for the partonic process du W^Z as a function of 
the angle 6 between the d quark and the Z boson in the du CM frame, at fixed energy 
-E'cM = 500 GeV and before any convolution with quark distribution functions. We have 
reported the different helicity contributions separately. As can be seen, the transverse 
component Utt shows the well-known radiation zero for cos^^ = (5'u,L+fi'd,L)/(5'u,L— fi'd.L) = 
— s^/(3c^) ~ —0.1 [ |2^, H, where g^^i^ and gd,L represent the Z-boson couplings to left- 
handed up and down quarks, respectively, and is strongly peaked in the forward and 
backward directions. The longitudinal contribution ctll is instead concentrated in the 
central region, at large angle of the Z-boson with respect to the incoming quarks. In- 
tegrating over the angle from 0° to 180°, one obtains the total cross sections shown in 
Fig. H (see the three curves on the left side) as a function of the energy. As expected, the 
dominant contribution is given by (Ttt, and above 300 GeV all polarized cross sections 
decrease with energy. 

The behaviour of the polarized cross sections depends, however, on the selected kine- 
matical region. If we consider the region of phase space characterized by a large transverse 
momentum of the Z boson, the relative size of the different helicity components and the 
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Figure 1: Lowest-order angular distributions for the process du W^Za' at Ecu = 
500 GeV. Here A, A' denote the transverse (T) or longitudinal (L) helicities. 




Figure 2: Born cross sections for the process du — > Za' as a function of Ecm with A, A' 
as in Fig. |l|. From left to right, the three legends refer to the left-side curves and to the 
right-side ones respectively, as explained in the text. 
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Figure 3: Lowest-order distributions in the invariants and yn defined in the text 
for the full process pp luil'V at y/s = 14TeV. Standard cuts and Pt{1'1') > 300 GeV 
are applied. 

shape of the curves change sensibly. As before, we plot the cross sections versus the CM 
energy but now for Pt(Z) > 300 GeV (see the three curves starting at around 600 GeV in 
Fig. D. In this case the LL contribution dominates at smaller energies, while the TT com- 
ponent increases with energy and takes over at high CM energies. This is due to the fact 
that the above-mentioned Pt cut translates into a minimum CM energy, -Ecm — 624 GeV, 
and limits the allowed range of the scattering angle of the Z boson. Hence, at low energies, 
the allowed angular region is strictly central and the LL component dominates. At larger 
energies, the allowed kinematical range increases by including smaller angles, and the TT 
contribution rapidly grows, soon overwhelming the LL part. 

Of course, one has to consider the additional effect due to the partonic distribution 
functions, which in turn decrease with increasing momentum fractions Xj and therefore 
with increasing CM energy \/S = y/xix^. The net result is shown in Fig. |^, where the 
distribution in the hard-scattering energy \/^ is plotted. Here and in the following we con- 
sider the full process pp ^ 4/, summed over all electron and muon final states, e^^i'le~^e~, 
e^'z/g/i+zi", /i^V^e+e~, and /i^ We have moreover applied our standard cuts as 

defined in Sect. ^ and the additional cut Pt{1'1') > 300 GeV on the reconstructed Z boson. 
As can be seen, despite the suppression resulting from the decrease of the parton distribu- 
tions with energy, roughly 50% of the contribution to the total cross section comes from 
the high-energy region a/5 > 1 TeV. We come back to this point later when discussing 
radiative corrections. 

As explained in Sect. |^, the DPA has proven to be a powerful tool for the computation 
of radiative corrections. In order to analyse, for the process considered here, the appli- 
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Figure 4: Bom cross section for the full process pp — luil'l' at ^/s = 14 TeV as a function 
of the cut on the transverse momentum of the reconstructed Z boson. Standard cuts are 
applied. 

cability of this approximation in a wide range of energies, we have plotted in Fig. § the 
tree- level cross section as a function of the Pt{1'1') cut. The first three curves represent, 
from top to bottom, the contribution of the pure doubly-resonant (DR) diagrams, the 
full result including all Feynman diagrams which contribute to the same final state (the 
number of diagrams is 10 in absence of identical particles in the final state, otherwise it 
doubles), and finally the DPA as defined in ( |3.1| ). 

If one looks at the difference between the complete result and the DPA, one can see 
that the discrepancy is rather remarkable. It amounts in fact to roughly 15% for Pt{1'1') 
cuts above 100 GeV. Let us note that, for this process, the commonly adopted narrow- 
width approximation, which corresponds to production x decay, differs from the DPA 
by less than 2% for the considered range of P|"*(/'F). So, depending on applied cuts 
and selected energy range, the narrow-width approximation can underestimate the exact 
result by roughly 13%. This points out the need of using the exact matrix element at 
lowest order, as specified in ( p.4|) . 

The second information one can extract from Fig. ^ is related to the contribution of 
non-DR diagrams. As shown by the dashed-line, the DR contribution (pp — >■ WZ — > 4/), 
which is lower than the exact result by about 1% around threshold, increases with energy 
relatively to the full result. For P™*(rr) = 300 GeV, the difference between the two 
cross sections is already of order 20%, and at very large energies the DR diagrams can 
even overestimate the result by a factor 2 or more. This effect is due to delicate gauge 
cancellations between DR and non-DR diagrams, which characterize the behaviour of off- 
shell cross sections in the high-energy regime. DR and non-DR diagrams do not constitute 
two separately gauge-independent subsets. Hence, the pure DR contribution cannot be 
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Figure 5: Bom cross section for the process pp — Ivil'V at y/s = 14TeV as a function 
of the upper cut M™* on the two invariant masses, M{ij), of the leptonic pairs which 
could reconstruct W and Z bosons, as explained in the text. Standard cuts and Pt{1'1') > 
800 GeV are applied. 

considered as a physical observable and the signal definition based on the diagrammatic 
approach and commonly adopted for example at LEP2 for WW and ZZ physics is not 
anymore adequate to describe di-boson production at the LHC in the high-Px region. 
The only sensible observable is the total contribution or the DPA which is a well-defined 
gauge-independent quantity. 

In order to investigate whether the difference between DR and full result is essentially 
due to the off-shellness of the gauge bosons as expected, we have then studied the effects 
of possible kinematical cuts. In Fig. ^, we have plotted the cross section for the extreme 
case Pt{1'1') = 800 GeV as a function of an upper cut M"^^^ applied on the two invariant 
masses M{ij) of the leptonic pairs which could reconstruct the Z and W bosons, M{ij) < 
Mv + M™*. We assume a lower cut M{ij) > My — 20 GeV, which is kept fixed in order 
to suppress the contribution from the virtual photon. As can be seen, for M™* = 20 GeV 
the difference between DR and the exact result reduces to the per-cent level. Also, both 
converge towards the DPA value, represented by the nearly flat dot-dashed line. It is 
quite obvious that cross sections computed in DPA are not sensible to this kind of cuts, 
as the gauge bosons are always considered on-shell except for the weakly cut-dependent 
factor [(p^ - M^)2 + F^M^]"^ - M|)^ + F|M|]"\ which reproduces the resonant 
peaking structure. 

Of course, a cut on the invariant mass of the lu pair is not physical since the longitudi- 
nal momentum of the neutrino is not directly measurable. We have therefore imposed the 



same kind of cuts, but using the transverse mass Mt(/z/) = yi?|(Zi/) — P|(Zi/) as physi- 
cal quantity instead of the M(/i/) invariant mass and releasing the lower cut on Mt^Iu). 
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The conclusion is similar. The DR contribution differs at the order of ten per cent from 
the previous case, as shown by the dotted curve in Fig. ^ The full calculation, which 
represents the true observable, and the DPA are instead rather insensitive to this change. 
In the following, we assume the additional kinematical cuts 

Mt(/z/) < Mw + 20 GeV, |M(/T') - Mz| < 20 GeV (4.1) 

under which exact result and DPA coincide at per-cent level, with a modest loss of signal, 
as shown in Fig. ^ where the lower solid line represents the full result after imposing the 
above-mentioned cuts. Since the exact cross section for the process pp 4/ is rather 
well approximated by the DPA if proper cuts are applied, we can safely adopt the DPA 
for computing electroweak radiative corrections. Let us notice, however, that electroweak 
radiative contributions are not much larger than 20% in the region of experimental sen- 
sitivity, as shown in the next section, and the DPA differs from the exact result by less 
than 15%. Therefore, without imposing the additional cuts ( ^.1| ), the error induced by 
use of the DPA in computing 0{a) contributions would give rise to an uncertainty of less 
than 3% on the total cross section, so well below the statistical accuracy. 



4.2 Effects of 0{cx) corrections 

In this subsection, we discuss the effect of leading-logarithmic electroweak virtual correc- 
tions to WZ production in DPA. First of all, one can see in Fig. ^ where the distributions 
in the reconstructed invariants and J\f\ are plotted, that the previously discussed 



conditions, under which the logarithmic high-energy approximation is valid, are well ful- 
filled for WZ production at high transverse momentum Pi[{ll). Both the hard-scattering 



invariant mass and y |f| = y\xiPi — where XjPj is the momentum of the parton 

from one of the protons (f corresponds to t or m depending on the partonic process), are 
in fact much larger than the boson masses. We have checked in addition that most part 
of the contribution to the cross section comes from the region where the scattering angle 
of the reconstructed Z boson in the WZ rest-frame is in the central range with respect to 
the beam. Finally, as to the ratios between the different invariants which appear in the 
logarithms, we have verified that the pure angular-dependent ones are mostly in the range 
1 < s/\r\ < 6, while s/M^ > 50 thus allowing to omit log^ (l^l/-?) type of logarithms up 
to an accuracy of a few per cent. 

As already mentioned in Sect. |^, we perform the computation of radiative corrections 
to the full process pp —>■ 4/ in DPA, using the complete expressions given in App. ^ i.e. 
implementing the full (SU(2)-transformed) Born matrix elements as in ( |A.12| ), (|A.13| ), 
( |A.18| ), and ([A.21|) -( [A.23|) . We have verified that the results obtained by making use 



of the high-energy approximation for the Born amplitudes given in ( |A.7| ) and ( |A.8| ) are 



in very good agreement. For all results given in the following, the difference between 
the two methods is ill fcLCt clt per-mille level. This comparison shows the reliability of 
the high-energy approximation for the Born matrix elements, under which the correction 
factor can be factorized and expressed in a very compact and simple form, leading to 
considerable decrease of CPU time. 

In order to discuss the basic structure of radiative corrections, we first consider the 
0{a) contributions to the partonic subprocess du — > W+Z. In Fig. ^ja we plot the 
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Figure 6: (a) Relative corrections All to the angular distributions for the process 
du WlZl- For each CM energy, the upper curves include the complete logarith- 
mic corrections, the lower curves only the angular-independent logarithms, (b) The same 
for the process du — > W^Zt. 



relative correction to the angular distribution of the longitudinal component All = 
(da^^/dcos^ - da^^j.Jdcos6)/{da^^j.Jdcos6) with a = a{dn W+Z) as a function of 
the angle 6 between the Z boson and the d quark in the du CM frame for the two energies 
Ecu = 0.5 TeV and 1.5 TeV. As can be seen, the LL part receives sizeable corrections, in 
particular in the central region 6 ~ 90° where ull is more enhanced. In order to pinpoint 
the effect of the angular-dependent contributions to the radiative corrections, in the same 
figure we have also plotted the two flat curves which include only angular-independent 
logarithms of s/M^. The difference between the two results for each CM energy shows 
the importance of taking into account leading and full subleading terms. There are in 
fact partial cancellations occurring between angular-dependent and angular-independent 
parts, which sizeably lower the overall corrections (see also Ref. [ p5[]). 

For the transverse part cttt the corresponding relative corrections Att are shown in 
Fig. ^ for two values of Ecu as in the previous figure. Here, radiative-correction effects 
are less pronounced compared to the LL case, especially at extreme angles where cttt 
receives its maximal contribution. The spikes in Fig. ^b originate from the radiation zero 
of the lowest-order cross section (see Fig. p, the absolute corrections behave smoothly 
everywhere. The angular behaviour of Att is more complex, compared with the longitu- 
dinal one. The dependence on the angle has in fact a two-fold origin. In addition to the 
angular-dependent logarithms [see (|A.18| ) and ( |A.19D ], there are angular- independent dou- 
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Figure 7: Relative correction to the total cross section for the full process pp — * luil'l' at 
^/s = 14 TeV as a function of the cut on the reconstructed Z-boson transverse momentum. 
Standard cuts are applied. 



ble logarithms log^(s/M^) with angular-dependent coefficients [see ( |A.13| ) and ( |A.15D ] 



which originate from the mixing of the final Z boson with the photon, induced by virtual 
soft-collinear W bosons. In Fig. |]b we have reported the total deviation Att, represented 
for each CM energy by the upper curves, and the partial contribution coming from the 
angular-independent logarithms of s/M^, given by the lower lines. As can be seen, also 
in this case, the correction factors proportional to the angular-dependent and angular- 
independent logarithms have opposite sign, leading as before to a reduction of the total 
correction. 

In order to show the effect of the electroweak radiative corrections on the complete 
process pp — > 4/, in Fig. 0we have plotted the 0{a) correction relative to the total Born 
cross section, A = (cr — (TBorn)/cBorn, as a function of the cut on the transverse momentum 
of the reconstructed Z boson, P^^*'{U). Our standard cuts are applied. As can be seen, 
the 0{a) contributions are negative and get larger with increasing P^^^, roughly going 
from —5% to —25% in the considered momentum range. The relatively large size of the 
radiative corrections, especially at energies which are at first sight rather modest (e.g. 
P™* > 250 GeV implies > 500 GeV), is mainly due to two combined effects. On one 
side, the longitudinal component of the cross section, ctll, which dominates at low values 
of the allowed energy range, as shown in Fig. ^ (right side), where the scattering angles 
are dominantly central, generates sizeable corrections. On the other hand, as can be seen 
in Fig. ^ the total cross section even for modest values of P™* receives a substantial 
contribution from the very-high-energy region, where cttt dominates. So, the generally 
smaller 0{a) contributions from the TT configuration, as compared with the LL ones, 
get enhanced by the higher values of the CM energy, and give globally additional sizeable 
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Figure 8: Rapidity distribution for the full process pp — > luil'l' at y/s = 14 TeV. Standard 
cuts and Pt^{VI') > 300 GeV are applied. The inset plot shows the difference between 0{a) 
and Born results normalized to the Born distribution. 



effects. As a consequence, the corrections to the total cross section are large because the 
Pt cut selects high-energy domains (since the only way to obtain a large Pt{U) is to have a 
large WZ invariant mass) and enhances the contributions coming from the central angular 
region. This shows that the size of the radiative contributions is strictly dependent on 
the applied cuts and the selected kinematical configurations. 

This is even more clearly visible in Fig. ||, where we have plotted the distribution 
in the difference between the rapidity of the reconstructed Z boson and of the charged 
lepton coming from the decay of the W, Ayzi = y{l'l') — y{i), at Born level (solid line) 
and including radiative corrections (dashed line). The rapidity is defined from the energy 
E and the longitudinal momentum Pi^hj y = 0.51og((i? + Pi^)/{E — Pl))- This variable, 
studied in Refs. [ ^ Q and defined in terms of direct observables, is symmetric around 
zero and shows a residual dip reflecting the approximate radiation zero of the angular 
distribution of the Born WZ production. The quantity Ayzi is in fact similar to the 
rapidity difference Ayzw = Vw — Vz considered in Ref . [ |[ , which is strictly related to 
the scattering angle, 6, of the Z boson in the WZ rest-frame. The definition of Ayzw 
and cos^ requires, however, the reconstruction of the unknown longitudinal momentum 
of the neutrino. Even if this can be derived by assuming the W boson to be on-shell [|26|, 



the two-fold ambiguity given by the two possible solutions for the neutrino longitudinal 
momentum spoils the radiation zero. Therefore, in order to extract informations about 
the angular dependence of the WZ process, it is preferable to use Ayzi [ 3 • 

The first information one can get from Fig. | is that the main contribution to the 
cross section originates from small values of Ayzi corresponding to central scattering 
angles. At the LHC, for the first time the statistics will be sufficient to experimentally 



15 



pp — s> luil'l' 


1 V / L J 


JjUill L J 


L J 


A [%] 


l/v/2L(TBorn \%] 


250 


1.716 


1.595 


-7.1 


5.4 


300 


0.899 


0.811 


-9.8 


7.5 


350 


0.503 


0.441 


-12.4 


10 


400 


0.296 


0.252 


-14.9 


13 


450 


0.181 


0.150 


-17.1 


16.6 


500 


0.114 


0.092 


-19.3 


20.9 



Table 1: Cross section for pp — > Ivil'l' for various values of P™^{1'1') 

test the behaviour due to the approximate radiation zero, which might be distorted by 
new-physics contributions. Figure ^ indicates that radiative effects are maximal at small 
rapidity separation, which is the region of stronger sensitivity to new physics. Moreover, 
owing to the applied cuts, these relatively large radiative contributions are not due to the 
suppression of the tree-level cross section and, being negative, they even slightly enhance 
the residual dip. 

The above-discussed effects should of course be compared with the expected experi- 
mental accuracy. In Table [I| we have listed the relative deviation A and the statistical 
error, estimated by assuming a luminosity L = 100 fb~^ for two experiments, for some 
P™* values. This comparison indicates that at high transverse momentum of the gauge 
bosons the virtual electroweak corrections are non-negligible and can be comparable with 
the experimental accuracy up to about 500 GeV. In this region the corrections range be- 
tween —5 and —20%. Whether or not they should be taken into account when performing 
analyses in this kinematical region depends of course on the available luminosity. Only 
in a high- luminosity run their effect will be relevant. 

5 W^7 production 

In this section, we extend our analysis to the process pp lui'j (/ = e,/i). This channel, 
proper for the measurement of the trilinear gauge-boson coupling WW7, can furnish 
complementary informations on the vertex structure of the SM when combined with the 
analysis of WZ production. As before, we consider the region of high CM energies of 
the hard scattering, where the sensitivity to new-physics effects is expected to be more 
enhanced, and the precise knowledge of the SM background can be then particularly 
useful. In the following we analyse the same set of variables used to discuss the WZ 
production process and the effect of the 0{a) electroweak corrections on them. 
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Figure 9: (a) Born cross section for the full process pp lui'-f at ^/s = 14TeV as a 
function of the cut on the photon transverse momentum. Standard cuts are applied, (b) 
Relative corrections to the angular distribution for the subprocess du — W^T- The two 
upper and lower curves for each energy are as in Fig. 

5.1 Born level 

The partonic subprocess du — >■ W"'"7 is dominated by the production of two transverse 
gauge bosons, whereas the remaining (LT) helicity configuration is suppressed by a factor 
Mw/.E'cM in the high-energy hmit. All features discussed in Sect. |]T]for the subprocess 
du W^Zt with transversally polarized gauge bosons qualitatively apply as well to W7 
production. The corresponding cross section is in fact strongly peaked in the forward and 
backward directions and presents a radiation zero for cos 6 = (Qu + Qd) / [Qn — Qd) = 1/3 
where 6 is the angle between the d quark and the photon. As to the general behaviour of 
the du —>■ W+7 process, we refer back to Figs. |I| and | and details given in the text. 

In spite of these similarities, W7 production presents, however, some different char- 
acteristics with respect to the WZ case. First of all, owing to the absence of any pure 
non-suppressed longitudinal components, there are no sizeable gauge cancellations in the 
total cross section of the full process pp — *• lui'y at high energy. The resonant contri- 
bution (pp W7 li^n) is always lower than the full result, also for high values of 
the cut on the photon transverse momentum Pt(7), and the difference between the two 
cross sections is below 3%. Therefore, one can still consider the pure resonant part as a 
useful definition of the W7 signal. Also, the single-pole approximation (SPA) defined in 
( p.2|) differs negligibly from the exact result, as shown in Fig. |a, where we have plotted 
the total cross section versus the cut applied on the photon transverse momentum. One 
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could directly use the SPA to compute radiative corrections at the per-cent level without 
imposing any additional cuts. However, for sake of uniformity, in the following we apply 
the same kind of cut Mt(/z/) < Mw + 20 GeV as used for the WZ process, under which 
resonant, full, and SPA cross sections converge to the same value within one per cent and 
with a negligible loss of signal. 



5.2 Effect of 0{a) corrections 

In this subsection, we study the effect of virtual electroweak radiative corrections on W7 
production in SPA. We consider first the partonic subprocess du W"''7. In Fig. 
we have plotted Att = {da'^'^ / d cos 9 — da^'^j.^/dcos9)/{da^Jj.^/dcos9) as a function of 
the angle 9 between the d quark and the photon. As can be seen, the behaviour of the 
0{a) contributions is quite similar to the WtZt case. Only the spikes, again due to 
the radiation zero, are correspondingly shifted, and reverse the shape of the curves with 
respect to the angle, compared with Fig. ^d. 

As in the previous case, the validity conditions of the high-energy logarithmic approx- 
imation for the radiative corrections are well satisfied by the complete process pp — /Z//7. 
The kinematical behaviour of W7 at high transverse momentum of the photon repro- 
duces in fact the same shape of the distributions as in Fig. |^. All invariants are then 
much larger than the boson masses, and at fixed P™*(7) the process receives considerable 
contributions from very high CM energies. We have moreover verified that, despite the 
radiation zero and the absence of any non-mass-suppressed longitudinal components, for 
large Pr^'^*(7) values (-P™*(7) ^ 250 GeV) most part of the contribution to the total cross 
section comes from the region of phase space where the photon is emitted at large angle 
with respect to the beam (see also Fig. pH]). 

In order to show the effect of radiative corrections on the full process pp li^n, we 
have plotted as before the 0{a) corrections relative to the Born cross section as a function 



of the cut on the transverse momentum of the photon in Fig. [T^. The overall behaviour 
is quite similar to the WZ case; the size of the radiative effects is, however, lower. 

This affects in the same way also the distribution in the difference between the rapidity 
of the photon and the charged lepton coming from W-boson decay [ Ayi^ = Vi — V-y, 
plotted in Fig. |1T[ Here, unlike in the previous case, the dip reflecting the radiation zero 
is much more pronounced, but the radiative corrections slightly decrease going towards 
Ayi^ = 0. Owing to the different location of the radiation zero, the radiative contribution 
for very small rapidity separation is still sizeable, but it does not get enhanced as for the 
WZ process, where in the same region the longitudinal component gives the dominant 
contribution. Hence for W7, despite the complex behaviour shown in Fig. ^d, which is 
merely due to the fictitious spikes, the radiative-correction effect is rather uniform in 
the angular range we consider, and leads to an overall rescaling of the Ayi^ distribution 
by roughly a factor 0.9. These effects could still mimic the behaviour of new-physics 
contributions. Their smaller size, compared with the WZ case, is compensated by the 
larger value of the overall cross section. Therefore, even if not extremely enhanced in the 
central rapidity range, radiative effects can become comparable with the statistical error. 

In Tab. 2 we compare the 0{a) relative correction A to the Born cross section with the 
expected experimental accuracy, assuming L = 100 fb~^ for two experiments, for different 
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Figure 10: Relative corrections to the cross section for the full process pp lui'-f at 
-y/s = 14 TeV as a function of the cut on the photon transverse momentum. Standard 
cuts are applied. 
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Figure 11: Rapidity distribution for the full process pp lui'f at \/s — 14 TeV. Standard 
cuts and Pt(7) > 300 GeV are apphed. The inset plot shows the difference between 0{a) 
and Born results normalized to the Born distribution. 
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values of P^^^{j). As one can see, radiative effects are very sensitive to P^^^{'y) and, despite 



pp lua 


Pt\i) [GeV] 


O-Born [fb] 


<7[fb] 


A [%] 


1/ ^J2La-Born [%] 


250 


5.810 


5.519 


-5.0 


2.9 


300 


3.180 


2.940 


-7.6 


4.0 


350 


1.832 


1.650 


-10.0 


5.2 


400 


1.100 


0.966 


-12.2 


6.7 


450 


0.684 


0.587 


-14.2 


8.6 


500 


0.437 


0.366 


-16.2 


10.7 


550 


0.285 


0.234 


-18.0 


13.2 


600 


0.190 


0.152 


-19.8 


16.2 


650 


0.129 


0.101 


-21.4 


19.6 


700 


0.089 


0.068 


-23.3 


23.7 



Table 2: Cross section for pp luij for various values of P™*(7). 

of the decrease of the cross section with increasing P™*(7), are larger than the statistical 
error for P^^^{'y) below 700 GeV, where they range from —5 to —23%. Moreover, they 
could be of some relevance also in a low-luminosity run (L = 30 fb"^) of the LHC, as they 
might become comparable with the experimental precision for P^'^*(7) < 400 GeV. 

6 Conclusion 

By means of a complete four-fermion calculation, we have examined WZ production in 
the purely leptonic channel at the LHC. An analogous computation has been performed 
for the W7 process followed by the leptonic W decay. We have given some examples of 
phenomenological analyses relevant to hadronic di-boson production in the high di-boson 
invariant-mass region. 

At tree level, we have found that, for processes involving WZ production, the dia- 
grammatic approach usually adopted to isolate the signal is not viable anymore at large 
transverse momentum of the reconstructed Z boson, owing to gauge cancellations. The 
doubly-resonant approximation can differ from the full result by tens of per cent in ex- 
perimentally relevant regions. The only sensible observable is the total contribution. 
Moreover the two commonly used approximations, narrow-width (i.e. production x decay) 
and leading pole approximation, can underestimate the exact result by about 10-15% at 
relatively modest energies, if no cuts are applied. 
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The primary aim of our analysis was to investigate the structure of virtual electroweak 
corrections and their effect on di-boson production processes at the LHC. The one-loop 
leading-logarithmic corrections to the full four-fermion or two-fermion-plus-photon pro- 
cess have been calculated in leading-pole approximation, neglecting non-factorizable cor- 
rections, and restricting oneself to the gauge-invariant leading-logarithmic corrections, 
which only contribute to the gauge-boson pair-production subprocess. We found that 
this approach constitutes a reliable approximation in the high-Px region at the LHC. 

In order to illustrate the behaviour and the size of 0{a) contributions, we have pre- 
sented different cross sections and distributions. In this study, we have not included the 
full QED radiative contributions, which involve also the emission of real photons and 
therefore depend on the detector resolution. We focused instead on the contributions of 
the leading electroweak logarithms resulting from above the electroweak scale. 

For WZ and W7 production processes, electroweak corrections turn out to be non- 
negligible in the high-energy region of the hard process, in particular for large transverse 
momentum and small rapidity separation of the reconstructed vector bosons, which is 
the kinematical range of maximal sensitivity to new-physics phenomena. Electroweak 
radiative effects lower the Born results by 5-20% in the region of experimental sensitivity. 
We have moreover shown that their size depends sensibly not only on the CM energy but 
also on the applied cuts and varies according to the selected observables and kinematical 
regions. Despite of the strong decrease of the cross section with increasing di-boson 
invariant mass, radiative effects can still be appreciable if compared with the expected 
experimental precision. This depends of course on the available luminosity. For WZ 
production, these effects are relevant for the high-luminosity run of the LHC. Owing to 
their larger overall cross section, W7 production processes can instead show a sensitivity 
to radiative effects also at low luminosity. 
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A Logarithmic electroweak corrections 

In this section, we present the analytical formulas for the logarithmic electroweak correc- 
tions to the polarized partonic subprocesses 

Mp-d) MPu) ^ W+^(pw) n^Apn), N = A,Z, (A.l) 



which can be derived from the general results given in Ref . [ |T3| . The label L indicates the 
left-handed chirality of the initial-state quarks (right-handed quarks are not considered 
since they cannot produce W bosons), and Aw.Af = 0, ±1 represent the gauge-boson 
helicities. The photon field is denoted by A in this appendix. The Mandelstam variables 
read s = {p^ + Pu)^, t = {Pd ~ Pw+)'^, and u = {Pd ~ PnY, where the momenta of the 
initial and final states are incoming and outgoing, respectively. In the high-energy limit. 
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we have i ~ — s(l + cos^)/2 and u ~ —5(1 — cos^)/2, where ^ is the angle between 
and pW, in the CM frame of the scattering quarks. 

For the calculation of the cross section ( |2.2| ) we need besides the process du — > W"'"A^ 
also its charge conjugate and the cross sections for exchanged initial quarks. These latter 
can be obtained from one another just by exchanging the invariants t ^ u. Owing to 
CP invariance, the charge-conjugate processes can be instead obtained from the initial 
processes by applying a parity transformation, 



M 



d(Pd)u(pfl) ^ W~(pw)iV(p^) = M d(pa)ufe) ^ W+(pw)iV(Piv) (A.2) 



with p = {E,—p) for p = {E,p). So, also in this case, the correction factors can be 
obtained from the same initial process (|A.1| ). The formulas we give in the following for 



the process ( |A.1| ) can therefore cover all contributions we need for the complete WZ and 
W7 production processes. 

The one-loop corrections are evaluated in the limit 

s ~ t ~ M > M^, (A.3) 

and we restrict ourselves to the combinations of gauge-boson helicities that are not mass- 
suppressed compared with in this limit. These correspond to the purely transverse and 
opposite final state (Aw, Aat) = (±, ^), which we denote by {Xw,^n) = (T, T), and, in 
the case of W^Z production, also to the purely longitudinal final state (Aw, Az) = (0, 0), 
which we denote by (Aw, Az) = (L, L). 

A.l One-loop corrections 

In the following, we present the results as relative corrections 

jr _ '^-^virt (Pd,Pu,PW,Pjv) .... 

to the Born matrix elements. A more detailed derivation can be found in Ref. [ p8[ . 



As shown in Ref. [ |13| , in the high-energy logarithmic approximation the longitudinal 
gauge bosons can be replaced by the corresponding would-be Goldstone bosons. There- 
fore, in our results for longitudinal final states (A = L), the substitutions Wl — > 0^ and 
Zl — ^ X have to be performed. 

The corrections (|A.4| ) are split as 

into leading (5^^^) and subleading (5^^^) contributions originating from soft-collinear 
gauge bosons, contributions 6'" that originate from collinear (or soft) gauge bosons and 
from wave-function renormalization, and contributions 6^^ that originate from parameter 
renormalization. All these corrections are evaluated in logarithmic approximation, i.e. 
including all terms that involve logarithms of the form log(s/M^) in the high-energy 
limit. More precisely, we restrict ourselves to 
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b) 





Figure 12: Dominant lowest-order diagrams for dLUL 



'X and dLUL — > W^A^t 



• the angular-independent double- and single-logarithmic corrections of the type 
alog^(s/M-^) and alog(s/M^), which involve only the ratio of the CM energy 
to the W-boson mass, 

• the double logarithms of the form a log (s/M^) log(M|/M^), and 

• the angular- dependent double logarithms of the type a log (s/M^) log (|r|/s), with 
f = t,u. 

For completeness we give also the analytic expressions of the double- and single- 
logarithmic corrections that contain logarithms log(M^/A^) and log(M^/mj), which 
involve the photon mass A or masses of light charged fermions.Q These contributions, 
denoted by L'^™ and Z'^™ in the following, are of pure electromagnetic origin and are not 
included into the numerical studies. 

The coefficients of the various logarithmic terms are expressed in terms of the eigen- 
values /^", or of the matrix components I^^,, of the generators^ 



Y 

-Q = — T^ 

^ 2 



■^W Y _|_ '^W rp3 



± iT2 
y2sw 



(A.6) 



where ^ = 1 - 4 = M^^/Ml 



A. 2 Born matrix elements in the high-energy Umit 



As input for the evaluation of the relative corrections (|A.4|) we need the Born matrix 
elements for the processes ( |A.1|) and the SU(2)-transformed Born amplitudes that we list 
in the following, restricting ourselves to the non- suppressed helicities. The corresponding 
amplitudes, are given in ( [A.7|) and ( |A.8| ) in high-energy approximation, i.e. omitting mass- 
suppressed terms. As we will see, this leads to very compact analytical expressions for 
the relative corrections ( |A.4| ). However, we recall that in the numerics instead of the 
high-energy approximations (|A.71 ) and ( |A.8|) the corresponding exact expressions have 



been used. As noted in Sect. 4.2, the difference is at the per-mille level. 



^This kind of contributions includes also energy-dependent double logarithms of the type 
alog(s/Af2^)log(M^/A2). 

detailed list of the gauge-group generators and of related quantities that are used in the following 
can be found in App. B of Ref. [ Oj. 
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For longitudinally polarized gauge bosons, we consider the Born matrix elements in- 
volving the corresponding would-be Goldstone bosons. These are dominated by the s- 
channel exchange of gauge bosons (see Fig. ^2^), and read 

-'^'Born ~ A) 2 £ ' 'Born ~ ^ QL xH 5 ' -'^'Born ~ 

-^Born - -e + gL-UL,dL. (A. ^ 

The production of transverse gauge bosons is dominated by the t- and u-channel contri- 
butions (see Fig. |l2|b) and gives 



. .gLgL^Af^-A^T _ ^2rAr'r7V / I , 4 . .gLgL-w+w^ _ e"^ At 

■^Born - e Iq^ ^ [j ^ ) ^^orn " ^^•^'> 

where f = t,u for g = d, u, respectively. In order to determine the relative corrections, the 
explicit dependence of the amplitudes Ag and A^, in ( [A.7| ) and ( |A.8| ), on the kinematics 
and on the helicities need not to be specified. 



A. 3 Leading soft— collinear corrections 

The angular-independent leading soft-collinear (LSC) corrections, which are given in 
Eqs. (3.6) and (3.7) of Ref. [ 0, depend on the eigenvalues 
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of the electroweak Casimir operator C"^™, on its components 
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in the neutral gauge-boson sector, as well as on the squared Z-boson couplings 
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For longitudinal and transverse gauge bosons we have 
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and 
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respectively, where the electromagnetic logarithms L'^™(s, A^, M^) are given by 
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For transverse final states, the non-diagonal components and C^za of electroweak 
Casimir operator require the evaluation of the transformed matrix elements with A^' ^ N. 
Using the high-energy approximation of the Born matrix elements ( [A.8| ), (|A.13|) can be 
written as 



a 



xLSC _ 



1 



log- 
log 



V3=dL,UL,W 



Ml 



2 2\ 



l/3=dL,UL,W" 

with the angular-dependent functions 
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A. 4 Subleading soft— collinear corrections 

The angular-dependent subleading soft-collinear (SSC) corrections are obtained 
by applying the formula (3.12) of Ref. [|TB[, to the crossing symmetric process 
^i.{Pd)'^i.{Pn)^x^{-P\N)Nxj,{-pN) 0, with ri2 = s, ri3 = t, and = u. This 
yields 
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In the cases A = L and A = T, the sums run over A^^ = XiH and Nl^^ = A, Z, respectively, 
and I^^ are defined in (B.23) and (B.27) of Ref . [ 13|. Using the SU(2)-transformed Born 
matrix elements given in ( |A.7D and (|A.8| ), we obtain 
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and is given in ( |A.16| ). 



A. 5 Single-logarithmic corrections 

The single-logarithmic corrections consist of the contributions 5^ and described in 
Sect. |A.1| . For longitudinally polarized final states, according to Eqs. (4.6) and (4.33) in 
Ref. [ the corrections 5^ read 
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and the parameter renormalization yields 
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where 



19/(6s^) is the one-loop coefficient of the SU(2) /3-function, and Aa{M, 



represents the running of the electromagnetic coupling constant from the scale to M^. 

If the final-state gauge bosons are transversely polarized then the log(s/M^) con- 
tributions in 5*" which are associated to the final gauge bosons cancel the log(s/M^) 
contributions originating from parameter renormalization, and according to Eqs. (4.6) 
and an analogue of (A. 11) in Ref. [ 1^ one obtains 
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where 6nz represents the Kronecker symbol. 
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